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Abstract 

Within the heavy quark expansion techniques for the heavy hadron weak decays we analytically 
compute the coefficient of the power suppressed dimension five chromo-magnetic operator at 
next-to-leading order of QCD perturbation theory with the full dependence on the final state 
quark mass. We present explicit expressions for the total width of inclusive semileptonic decays 
including the power suppressed terms and for a few moments of decay differential distributions. 
One of the important phenomenological applications of our results is precision analysis of the 
decays of bottom mesons to charmed final states and extraction of the numerical value for the 
CKM matrix entry \V c b\. 



1 Introduction 


Presently the Standard Model of fundamental interactions is being thoroughly tested experimentally 
at colliders, but no definite signs of New Physics have been detected beyond the framework of the 
Standard Model. Neither new particles have been explicitly seen nor any significant deviations from 
the Standard Model values in the loop sensitive Wilson coefficients for flavor changing observables 
have been determined in high precision data (as a review, see e.g. [T]). Thus, the Standard Model 
has successfully passed all tests in the areas where it is certainly valid as a low-energy effective 
theory. 

However, there is definitely life beyond SM. Some new phenomena - like neutrino masses and 
mixing - can be readily incorporated in a rather straightforward manner to extensions of SM. 
The other new effects - like dark matter - are of cosmological nature and related to still poorly 
understood realm of gravity and, strictly speaking, are outside the physics of the standard model 
domain. Nevertheless it seems certain that the scale of the traditionally expected extensions of the 
standard model - like supersymmetry or extra dimensions - has definitely moved from few TeV 
region to a higher one in energies that can make it unreachable at accelerators in foreseeing future, 
e.g. [2]. Since the New Physics scale moved higher the direct observation of new physics phenomena 
will not probably be explicit even at new machines (still one should wait for the results of the 14 TeV 
run of LHC!). In case that nothing will be seen the new phenomena beyond the standard model 
(if any at all!) can only be identified through detecting slight discrepancies between theoretical 
predictions within the SM and precision measurements at low energy with available tools. 

Accurate theoretical predictions within the SM are of crucial importance in such a scenario. For 
these predictions to be reliable one first needs the precise numerical values for the key parameters of 
the SM itself. The least precisely quantitatively known sector of the SM is a quark flavor one where 
the quark Yukawa couplings to the Higgs held are not well known numerically. In the standard model 
they translate into the mixing angles between generations gathered in the CKM matrix and the 
vacuum expectation value of the Higgs held. The latter can be determined from the leptonic sector. 
Note that the havor sector is also a most promising place in investigating the Higgs mechanism 
that is definitely of an effective origin and probably will be modified in future as the presence of 
a fundamental scalar in the “final” theory does not look convincing. All in all the havor physics 
of quarks is the promising place to search for new physics and should be thoroughly studied (see, 
e.g. PI gj). 

While the quark weak decays are mediated by the charged weak currents at tree level, which 
are believed not to have sizable contributions of possible new physics, their study is of importance 
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for precise determination of the numerical values of the CKM matrix elements. However, obtaining 
solid theoretical predictions for processes with quarks at the fundamental level requires the use of 
genuinely nonperturbative computational methods like QCD lattice calculations since eventually 
one has to make prediction for the experimental quantities that include hadrons and cannot be 
described in perturbation theory of QCD due to confinement. This is principal part of the problem 
but there is also a pure technical part. Even if the direct computation in terms of quarks would be 
relevant to the world of hadrons that partly can be made possible by choosing proper observables one 
will still face the problem of computational complexity of the calculation with sufficient accuracy 
that requires a rather large order of perturbation theory. The example is the description of the 
process b —» sy. 

Taking just the parton level of computation for hadronic processes one makes the technical 
part equivalent to that of the leptonic calculations where the benchmark level for the technical 
part of the computation is the evaluation of the muon lifetime. The muon decay is a source for 
the determination of the Fermi constant Gf with high accuracy from a leptonic sector. First 
radiative corrections have been computed long time ago [3 EJ. To match the precision of the 
present experimental data for muon lifetime, the theoretical calculations have to be performed with 
very high accuracy. In this case the calculations are feasible, since the purely leptonic decays are 
well described within perturbation theory and the expansion parameter a ~ 1/137 is small. The 
latest theoretical result includes the second order (NNLO) radiative corrections in the fine structure 
constant expansion [7j 
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Here is the muon mass. The numerical value for the electron mass m e is set to zero everywhere 
but in the expression for the expansion parameter 
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The expressions with account for nonvanishing electron mass are known. The quantity /\r had = 
—0.042 ± 0.002 is the hadronic contribution that is known with uncertainty of about 5%. It cannot 
be computed from first principles for light quarks and is obtained by integrating the experimental 
data for the photon vacuum polarization. Note that the similar situation emerges with precision 
analysis one of the key leptonic observable - the muon anomalous magnetic moment g — 2. At 
present the hadronic contributions related to light quarks give the main uncertainty of theoretical 
prediction (e.g. my It is a general feature that quark sector influences even pure leptonic processes 
if the required accuracy is high enough, e.g. [10]. 
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Eq. (JTj) results in an C?(lppm) accuracy of the theoretical expression for the lifetime that is 
competitive for precision comparison with modern experimental data. As for the quark sector is 
concerned there is a good set of data for s —> u weak transitions that corresponds to K —> 7 \ev e 
decays at the hadron level, but it is hopeless to compute the related rate theoretically at present 
because of strong infrared problems in theoretical treatment of reactions with light hadrons. 

For heavy hadrons the theoretical treatment of the decays is however possible because the large 
mass of the heavy quark constitutes a perturbative scale that is much larger than A qcd- The leading 
logarithmic effects related to that scale have been discussed long ago m • Later there have been 
created a framework for the possibility for an expansion in powers of A qcd/^q where ttiq is the 
quark mass and A qcd ~ 500 MeV is a typical hadronic scale nsnsuii]. Top quarks do not form 
mesons due to their short lifetime, charmed mesons are probably not heavy enough, rendering the 
convergence in the inverse mass marginal, but the case of bottom-meson decays is certainly tractable 
in this way and thus has been intensively studied. The technique is applicable to b —>• u and b —> c 
transition and both to scmileptonic and purely hadronic inclusive decays. For definiteness, we will 
stick to semileptonic b —> c decays. 

In the present paper we analytically compute the coefficient of a power suppressed dimension 
five chromo-magnetic operator at next-to-leading order of QCD perturbation theory with the full 
dependence on the final state quark mass. The results of the analogous computation in the mass¬ 
less limit for the final state quark have been presented earlier in ref. P3j Here we present explicit 
expressions for the total width of inclusive semileptonic decays and few moments of differential 
distributions with full dependence on the final state quark mass. One of the important phenomeno¬ 
logical applications of our results is precision analysis of the decays of bottom mesons to charmed 
final states and an extraction of the numerical value for the CKM matrix entry V c f ,[. 

The paper is organized as follows. In the next section we give a general representation for the 
decay width of a heavy hadron in a form suitable for computation in QCD. In Sect. we give 
necessary basics of Heavy Quark Effective Theory (HQET) that is a working tool for the present 
calculation. In Sect. [4] we write down the Heavy Quark Expansion (HQE) for the decay rate. 
The actual computation and results are described in Sect. [5j In Appendices we give the explicit 
expressions for our master integrals and some long analytical expressions for the coefficients of HQE. 
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2 QCD representation for the decay rate 


It is difficult to compute an hadronic decay rate since the underlying theory of strong interactions - 
QCD - is formulated in terms of quarks and the hadrons only appear in the strong coupling regime 
as bound states. Therefore one can use either numerical calculation on the lattice or find special 
observables for which perturbation theory calculation is feasible in some form. Such observables are 
inclusive ones since the sum over hadronic states can be related to the sum over the quark-gluon 
states using unitarity of the theory. In case the initial state is treatable in perturbation theory, 
i.e. it is a leptonic one as in e + e _ -annihilation into hadrons or hadronic r-lepton decays then the 
results can be uniquely obtained in perturbation theory. In cases when the initial state is hadronic, 
i.e. it is non-treatable in perturbation theory, one uses a factorization idea - to separate scales 
and compute the short distance effects in perturbation theory while long distance properties are 
coded in hadronic matrix elements. The famous example of the latter approach is the analysis of 
deep inelastic scattering of leptons on hadrons. The analogue of deep inelastic scattering in heavy 
quark physics is inclusive decays of heavy hadrons. One can use either fully hadronic (non-leptonic) 
decays or semileptonic ones. The number of experimental observables in inclusive hadronic decays 
is however limited to basically the total rate of the process. In semileptonic decays the presence of 
leptons in the final states gives more kinematical flexibility still retaining the rigorous theoretical 
description of the process. 

The low-energy effective Lagrangian £ e g for semileptonic b —> c transitions is a Fermi four- 
fermion one 

C c s = 2v / 2G_pV r C 6(6L7/i c n)(^7 /i ^n) + h- c - (2) 

with left-handed fermion fields. The numerical value for Fermi constant Gf is determined from pure 
leptonic weak processes and known with high precision. The mixing angle V c b is the main interest 
in decay measurements with hadronic initial states [16]. The precision analysis of such processes is 
important both for the flavor sector and Higgs mechanism investigations in search for new physics. 

Using unitarity of the 5-matrix the inclusive decay rate B —>■ is obtained from taking the 

absorptive part of the forward matrix element of the transition operator T ra that is the second 
order term of the perturbation theory expansion in the interaction Lagrangian £ c q, 

T = ij dxT {C cS (x)C eS ( 0)} . (3) 

Note that the transition operator T is a non-local functional of the particle fields and is given 
by the integral over all possible scales. There is no much hope to handle such an operator in 
QCD that includes all scales as well and no large parameter is available in case of the two-point 
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correlator in Eq. (J3|. However, one can hope that some transitions or matrix elements are still 
short distance dominated even if it is not a universal feature of the correlator given in Eq. (|3]) itself 
and may depend on external states. For light hadrons (like kaon) it is definitely not the case and 
taking matrix elements cannot help in isolating a short distance dominant part of the correlator 
in ([3]) . Heavy hadrons have an additional simplification that makes the computation of some matrix 
elements possible in perturbation theory by separating the scales involved. 

The idea is that when taking a matrix clement over a heavy hadron containing a heavy quark 
with mass m,Q A qcd the correlator does acquire a large internal scale, mq , that enables scale 
separation. For actual separation of scales one applies the operator product expansion (OPE) tech¬ 
niques. These ideas are formalized through the notion of effective theories. Within the heavy hadron 
with momentum pu and mass Mh the large part of the momentum is due to a pure kinematical 
contribution of the heavy quark pu = rriQV + A with v = Ph/Mh being the velocity of the hadron 
and A is related to the light degrees of freedom and interactions between them and the heavy quark. 
One can already extract the factor related to the large quark part of the momentum explicitly at 
the level of field variables when afterwards the matrix element over a heavy hadron is taken. The 
heavy quark field can be separated into the fast oscillating phase and a slow changing field h v (x) 
with a typical momentum of order A ~ J^qcd 

Q(x) ~ e-< m ^ x h v (x) . (4) 


The velocity v = Ph/M^ is finite in the limit of infinitely heavy quarks itiq A qcd- This program 
is realized within the effective theory for heavy quarks. In order to make the dependence of the 
decay width on the heavy quark mass rriQ explicit and to build up an expansion in A qcd/^Qi one 
matches a time-ordered product of full QCD operators in entering to the transition operator T onto 
an expansion in terms of Heavy Quark Effective Theory (HQET) [T8j 29] • Presently the Heavy 
Quark Expansion in inclusive semileptonic b —> c transitions provides a level of theoretical precision 
in the prediction of the total inclusive rate for B —> X c tvn within two percent. The structure of the 
HQE is given by [20] 


T(B X c iv t ) = r°\V cb \ 
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where T 0 = G 2 F m\/{ 1927T 3 ) and mb is the Aquark mass. The precise definition and the proper 
choice of the most suitable mass parameter for the heavy quark field is extensively discussed in 
the literature. The power suppressed terms are given by the forward matrix elements of the local 
operators of growing dimensionality in HQET over the heavy hadron state. Their numerical values 
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are determined by the corresponding power of the QCD infrared parameter of A qcd- These are 
nonperturbative quantities either to be computed within some non-perturbative techniques such 
as lattice QCD or to be fitted to experimental data. The kinetic energy parameter is given 
by the nonrelativistic kinetic energy operator of the heavy quark within the heavy hadron. The 
chromo-magnetic parameter /Iq is given by the matrix element of the magnetic dipole operator. 
These two operators give the leading power suppressed contribution and were intensively studied. 
The higher order power suppressed terms are becoming important at present as the experimental 
data improves. The parameter p 3 describes the contribution of dimension six operators that are 
Darwin term and spin-orbit interaction. The general parameter p 4 is a contribution of a rather 
large number of dimension seven operators [21] . The coefficients a* are functions of the quark and 
lepton masses and have a perturbative expansion in the strong coupling constant a s (mb)- The 
leading term coefficient ao is known analytically to O (a^) precision in the massless limit of the 
final state quark m At this order the mass corrections have been analytically accounted for the 
total width as an expansion in final fermion mass in ref. [23j and for the differential distribution 
numerically in [23j. The coefficient of the kinetic energy parameter is linked to the coefficient 
ao by Lorentz invariance, see the explicit analysis in [23], The NLO correction to the coefficient 
of the chromo-magnetic parameter 02 has been investigated recently in [26j where the differential 
distribution has been computed and the total decay rate has been then obtained by a process of 
numerical integration over the phase space. The a s correction to the chromo-magnetic parameter 
coefficient 02 has been analytically computed in ref. p~5j in the massless limit. Here we give the 
result with full mass dependence in analytical form. Our calculation of the coefficient 02 is in fact 
a matching computation between QCD and HQET. For this reason we present some facts about 
HQET relevant for our discussion in the next section. 

3 Basics of HQET 

A heavy quark near its mass-shell is described by a field h v (x) which is a remnant of the whole QCD 
fermion held Q(x). In fact, it effectively contains only large components of the Dirac bi-spinor that 
describe the quark and not the antiquark. One achieves the separation of the components by using 
the projector P + = (l+^)/2 where v is the external velocity that determines the remnant fields h v (x) 
and the whole construction of HQET. Note that obtaining HQET as the effective theory from QCD 
is very close in spirit to the well known procedure of obtaining the nonrelativistic limit of QCD or, 
earlier, QED. The held variables and Lagrangians are just the same in both nonrelativistic QCD and 
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HQET. The quark velocity v is fixed in the presence of the heavy hadron by its momentum. Usually 
the common choice for the velocity is v = Pr/Mr. The behavior of time and space components of 
the formal Lorentz four-tensors differs in HQET. It is useful to split a four-vector in longitudinal 
and transverse parts, namely = v ,J {vp) +p^_. The covariant derivative of QCD is 7 r M = id^ + 
with the splitting 7 = v^{im) + 7Tj(. 

The quantity h v is the heavy-quark field entering the HQET Lagrangian mm- The effective 
Lagrangian of HQET can be obtained in a concise form at tree level by integrating out the P_ part 
of the heavy quark field H v , H v = P-H v , with the result 


£-hqet — h v (7iv)h v + h v f A 


hy . 


(5) 


■ 2m + 7 tv ' 

Here the first term is just the residual energy of the quark while the second one describes the 
effects of the removed (integrated out) antiquark. It is non-local that is the price for integrating 
the antiquark out. In the limit m nv one can expand the second term in a series in the inverse 
large mass and obtain a local Lagrangian up to a given order in the mass expansion 


- - 1 - - 71V 
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( 6 ) 


It is inconvenient to have time derivatives in a term that is formally a correction since then the 
fields h v are not correctly canonically normalized. Therefore the redefinition of the fields is used to 
remove time derivatives 


h v —y ( 1 + 
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and get the Lagrangian for the new modes h v (for which we retain the same notation though) in 
the form 


£hqet ~ ° v + Tn^°' n 
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with 

C mSLg (p) = 1 H —Icp + Ca + hi T -^ (9) 

being the coefficient of chromo-magnetic operator O g including the QCD radiative correction of 
the order a s cat For new modes h v the terms of the order 0(l/ml) in the Lagrangian contain no 
time derivative [I9l [28] . Here we introduced the notation used below. The quantity O v = h v vnh v is 
the leading power energy operator that is independent of the heavy quark mass and spin and gives 
the famous spin-flavor symmetry of HQET. The quantity O n = h v ii\h v is a kinetic energy operator 
and O g = h v a flp [R^_,TT h []/2h v is a chromo-magnetic operator. They constitute classical subleading 
power operators. Higher terms are given by the operator 0 3 = h v [f ± , [f ± ,7rv]]h v that can further 
be converted into a linear combination of the Darwin Or = h v [ [7r((, nv]\h v and spin-orbit term 








Osl = [ 7r ±; 7rw ]]^! ^3 — cdOd + cslOsl with coefficients cd,csl known at the next-to- 

leading order of perturbative expansion in the strong coupling constant. The discussion of order 
1 /m 2 Q terms in the HQET Lagrangian is relevant for our computation because of the necessity to 
precisely fix the definition of the fields h v entering the heavy quark expansion. 


4 HQE for the width correlator 

For further convenience we introduce a normalized transition operator T through the relation 

imT = r°|F c 6 | 2 f. (io) 

With the use of heavy quark effective theory the heavy quark expansion is simply a matching from 
QCD to HQET 


O v 0 7T 

T=C 0 O 0 + C v ^ + C n -^ 
mb 2mf 


+ a 
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G 
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2 ml 


( 11 ) 


The local operators Oi in the expansion (11) are ordered by their dimensionality O o = h v h v , 


O v = h v virh v , O n = hyTT^hy, Oq = ]h v . The coefficients of these operators are obtained 

by matching the relevant matrix elements between QCD and HQET. Note that after taking a matrix 
element over the hadronic state (like the H-meson) one can use equations of motion for HQET fields 
h v to eliminate the operator O v . By the same token there is an operator 0 5 = h v (v7i) 2 h v that is of 
higher order in the large mass expansion after going on shell using equations of motion of HQET. 
Thus, the expansion ( [TT| is a matching relation from QCD to HQET with proper operators up 
to dimension five with the corresponding coefficient functions. The coefficients are independent of 
external states and one can take them at will. We take a heavy quark on shell and gluons as external 
states for matching to QCD. 

Note that one can use the full QCD fields for the heavy quark expansion expansion as well. 
However the choice of the proper basis of operators is not so straightforward as in HQET. Still it is 
convenient to choose the local operator bfb defined in full QCD as a leading term of heavy quark 
expansion [29]. Indeed, the current try^b is conserved and its forward matrix element with hadronic 
states is absolutely normalized. For implementing this setup one needs an expansion (matching) of 
a full QCD local operator b'fh in HQE through HQET operators. The expansion reads 


bjb = Oq-C^ 


2 ml 


+ ^ G 2mf + 0(^QCD/ m b) 


( 12 ) 


up to necessary order in the strong coupling a s . The coefficient of the leading power operator O 0 
has no radiative corrections and the kinetic operator has the coefficient related to the leading one 
due to Lorentz (reparameterization) invariance. 
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Substituting the expansion (12) into eq. pT| ) one obtains after using the equation of motion for 
the operator O v in the forward matrix elements 


f = o> 1 b/b - ^|} + {-c'.cW/i) + c G - c G c 0 } (13) 

Note that for phenomenological applications the numerical value for the chromo-magnetic moment 
parameter fi 2 ,, related to the forward matrix element of the operator O g , is usually taken from the 
mass splitting between the pseudoscalar and vector ground-state mesons. The mass difference of 
bottom mesons m 2 B * — m 2 B = A m 2 B = 0.49 GeV 2 is given by 


2 M, 


B 


-C ra&g {n){B(pB)\0 G (n)\B(p B )) = - Am 


B 


(14) 


(up to higher order 1/mg corrections) where we use the relativistic normalization of states. There¬ 
fore the coefficient in front of the renormalization group invariant combination C majg (p)0 G (/j,) can 
be useful. In such normalization one gets after taking the forward matrix element of the expansion 


in Eq. (13) the representation 


T{B X c v t i) = r°|Ecfc| 2 Co 1 + 


2 ml 
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(15) 


5 Description of the calculation and results 


5.1 Generalities and techniques 

The matching procedure consists in computing matrix elements with partonic states (on-shell quarks 
and gluons) at both sides of the expansion ( [II] ). The coefficient function Cq of the dimension three 
operator h v h v determines the total width of the heavy quark and at the same time the leading 
contribution to the width of a bottom hadron with HQE technique. At NLO the calculation of the 
transition operator T in ([ 3 J requires to consider three-loop diagrams with external heavy quark 
lines on shell. The leading order result is well known and requires the calculation of the two-loop 
Feynman integrals of the simplest topology - the sunset type ones |fJD]. At the NLO level one needs 
the on-shell three-loop integrals with massive lines due to the massive c-quark. In Fig. [I] we show 
some typical three-loop diagrams both for the partonic part and power corrections of the decay 
rate. 

The computation has been performed in dimensional regularization used for both ultraviolet 
and infrared singularities. We used the systems of symbolic manipulations REDUCE ra and 
Mathematica | 32 | with original codes written for the calculation. The package FeynCalc [33] is used 
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Figure 1: Perturbation theory diagrams for the matching computation at NLO level, (left) - 
partonic type, right - power correction type: insertion of an external gluon 

for manipulating Dirac matrices and four vectors under Mathematica. The reduction to master 
integrals has been done within the integration by parts technique [33]. The original codes have 
been used for most of the diagrams and then the program LiteRed [33] has been used for checking 
and further application to complicated diagrams. The master integrals have been computed directly 
and then checked with the program HypExp [36] • The renormalization is performed on-shell by the 
multiplication of the bare (direct from diagrams) results by the on-shell renormalization constant 

^ = 1 - c 'Tr(7 + 31 n © +4 )' (16) 

It is convenient to fix the normalization point to the 6-quark mass /j = mb in the practical compu¬ 
tation. The /i-dependence can be easily restored from the knowledge of anomalous dimensions. 

We present and discuss the obtained results below. 

5.2 The leading power coefficient C 0 : partonic width 

By using the described methods we reproduce the known result for the heavy quark width which is 
given by the contribution of the leading operator Oq. The coefficient Co is 

C 0 = Cq° + Cf—Cq LO (17) 

7T 

where the LO contribution reads 

Cq° = 1 — 8r — 12r 2 ln(r) + 8r 3 — r 4 (18) 
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and the NLO contribution reads 


rtNLO 
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(19) 


+ ^—10 — 45r + - r 2 —--r 3 J r ln(r) + ^—18 —— J r 2 ln 2 (r) 

+ (2 + 60r 2 + 2r 4 ) ln(l — r) ln(r) + (l + 16r 2 + r 4 ) (3Li 2 (r) — 7t 2 /2) 

+16r 3/2 (l + r) ^7r 2 - 4 (Li 2 (^/r) — Li 2 (—s/r)) + 21n(r) In j + 

with r = m 2 c /m 2 . Here Li 2 (r) is polylogarithm, Li 2 (r) = ^2 n r n /n 2 . The combination 

r 1/2 (ir 2 - 4 (Li 2 (v / r) - Li 2 (-\/r)} + 21n(r) In | + (20) 

is a part of one master integral in the computation and it always appears in thios form. It contains 
a specific odd contribution r 1 / 2 7r 2 while the rest is in fact formally even in m c . The analytical 


expression at NLO in Eq. (19) has been first given by Nir 

The behavior near the border of the decay phase space (r ~ 1) of the NLO correction 


riNLO 

L o 


(r^l) = -(l-r) 5 + 0((l-rf) 


( 21 ) 


is similar to that of the LO which is 


C 0 IO (r^l) = ?(l-,r) 5 + O((l-r)«). (22) 

A typical feature of the result at next-to-leading order is the presence of odd powers of the charm 
quark mass like r 3 / 2 . Of course, it does not mean that there is a symmetry m c —> —m c . At the 
small mass limit r —> 0 only the simplest term of such structure 7r 2 r 3//2 survives with a rather large 
coefficient. 

We define the bottom quark mass being a pole one because it is convenient for computing the 
relevant matrix elements in QCD with on-shell quark states. The definition of charmed quark mass 
can be either the pole scheme or MS-scheme one. The relation between the two definitions up to 
necessary order is 

K° l = { 1 + C ^ ( 3 ln ^ + 4 ) ) • ( 23 ) 

The numerical value for the charmed quark mass is best known in the MS-scheme [38l [39] . It is 
rather small and cannot be be perturbatively cast into the pole mass scheme with any reliable 
control over uncertainties due to convergence of perturbation series expansion ffl. The numerical 
value for the bottom quark mass has been discussed in the literature for a long time and many 
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Figure 2: The mass dependence of the coefficient Co(r) in a pole mass scheme for m c : solid line - 
LO, dashed-dotted line - NLO normalized to Cq LO ( 0): C^ LO (r) = C^ LO (r)/C^ LO (0), Co°(r) = 
C 0 LO (r)/C 0 LO ( 0) . 

estimates are available. Also there is a extensive discussion which particular scheme of defining the 
quark mass parameter which is the most suitable for this particular observable m 133- 

In Fig. [ 2 ] we give the plot of the coefficient Cq 0 [t ) and also the normalized next-to-leading 
coefficient C,f i0 (r) in the pole mass scheme for m c . 

In Fig. [ 3 ] we give the plot of the mass dependence of the coefficient Cq LO {t ) in different mass 
schemes for m c . 

In the small mass limit for the charmed quark one finds 



We have computed the results for the coefficient Co in massless limit, Co(0), independently that 
serves partly as a check of our full mass calculation. 

The relative magnitude of the NLO contribution at a typical value of mass ratio r = 0.07 is 

ry ry 

C 0 (0.07) = 0.60 - 0.78Cf— = 0.6(1 - C F — 1.31) (25) 

7 r 7T 

while in massless limit it is 

ry 

C 0 (0) = 1 — C F — 1.8. (26) 

7T 

The numerical value for the bottom quark mass m 5 is important for phenomenological applications 
and discussed in the literature (see, e.g. [II]). The dependence on the charm quark mass is essential 
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Figure 3: Mass dependence of the coefficient Cq LO {t) in pole and MS schemes with /j = rn b and 
li = m c . 

but still it follows mainly the pattern of that at leading order. This similarity supports the idea 
of ref. p2] that the computation in massless limit can be useful for physical applications as the 
normalization and the extrapolation with the leading order massive result can be a reasonable 
approximation for the mass dependence at NLO. We will see how it works or does not work for 
other coefficients later. 

5.3 The vD -operator coefficient C v 

Here we present the result for the coefficient C v which is an auxiliary quantity in our approach since 
the operator is reexpressed through the other contributions at the level of matrix elements. The 
coefficient C v is singled out by taking the matrix element between 6-quarks on shell and one gluon 
with vanishing momentum and longitudinal polarization, i.e. the gluon held is chosen on the form 
Afj, = v tl {vA). Here A p is a matrix in color space A p = A a ^t a . The result for the coefficient C v 

C v = C^° + C f —C; nlo (27) 

7T 

reads 


C™ = 5 — 24r — 12r 2 ln(r) + 24r 2 — 8r 3 + 3r 4 . 
In Fig. [4] we plot the charmed quark mass dependence of C v . 
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Figure 4: Mass dependence of the coefficient C v (r) 


For the NLO part C^ LO we give explicitly only the expression for the small r expansion. The 
structure of the whole contribution is very similar to that of Cq L °. The expression is rather long 
and given in Appendix [B] The small mass expansion reads 


C, 


NLO 


^ + 48r — 8n 2 r 3 / 2 + O ( V 2 In 2 r) . 


(28) 


The leading term of the expression coincides with the independent computation in the massless 
limit done in [15] 

a,| r=0 = 5 + CA^{-— -y} . (29) 

As for the mass dependence of the coefficient C v . for the typical value of r = 0.07 one finds 


O' / ry 

67(0.07) = 3.6 - 3.8 C F — = 3.6 (1 - C F — 1.1 

7T V 7T 


while in the massless limit one has 


Of e 


0,(0) = 5(1 -0—1.2) 


7T 


(30) 


(31) 


One sees again a rather reasonable accuracy for the mass dependence extrapolation at NLO. 

The coefficient C v has no Ca color structure, it contains only the CV Casimir invariant. This 
property matches the possibility to compute this coefficient using a small momentum expansion near 
the quark mass shell, p = mv + k. Still, an explicit cancellation of the contribution proportional to 
the color structure Ca and cancellation of poles with the same renormalization constant Z® 5 shown 


in Eq. (16) is a powerful check of the final result. 
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The large m c behavior at the border of phase space is 


Cf LO (r -)■ 1) = -3(1 - r) 4 + 0((1 - r) 5 ) (32) 

and 

C/ LO (r ->• 1) = 4(1 - r) 4 + 0(( 1 - r) 5 ). (33) 


5.4 The coefficient Cq — CqCq = Cq: 
chromo-magnetic operator 


For the chromo-magnetic operator coefficient we directly compute the difference between contribu¬ 


tions to the width correlator in Eq. (11) and the local bfb operator in Eq. (12) multiplied by the 
leading power coefficient Co(r), Cq = Cq — CqCq . We write this coefficient as leading order term 
and radiative correction in the form 


Cq{t) 


Qyr,LO 


(r) + 


a* 


7T 


^yC (J r ’NLO,A 


(r) + C F C\ 


.r,NLO 



(34) 


where the NLO coefficient is separated into two color structures with Ca and C F color group 
invariants. In Fig. [5] we present the plot of the mass dependence for the coefficient of the chromo- 
magnetic operator for QCD with Ca = 3 and Cj? = 4/3. One sees that the mass dependence of 



r = rn 2 c l m 2 b 


Figure 5: Mass dependence of the coefficient Cq = Cq — C 0 Cq with // = m b 
Cq at NLO is much sharper than in previous cases. This is unexpected and makes the conjecture 
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about a uniform phase space suppression for the coefficients less accurate. The explicit leading 
order expression reads 


C r a LO = 2-16r - 24r 2 ln(r) + 16 r 6 - 2r 4 = 2 C { 


LO 


(35) 


The NLO coefficients with full mass dependence are too long, whereas the expanded results are 


n r,NLO,A 


87 T 2 y/r ( 9/ 2:/r 2 \ 7T 2 49 

-h r ( ln 2 (r) — 25 ln(r) H-25 )- 1 -, 

3 V 3 ) 918’ 

n r,NLO,F 327 T 2 y/r , 2 . 2 , Q1 \ ^7r 2 47 

Q = - 3 -hr (-4In (r) + 681n(r) - 4 tt +21)---— . 


(36) 


At the border of phase space we obtain 




r,NLO 

G 


(r -H) = CV(1 - r) 4 + - (2 C f - 3C a ] (1 - rf + 0(( 1 - r) e ) 


(37) 


and 


Q LO (r->l) = -(l-r) 5 + 0((l-r) 6 ). 


In the massless limit the Cq coefficient is given by 




(38) 


(39) 


This result has been independently determined by the direct computation using the technology 
developed for the massless case. 


5.5 Coefficient = — C v + C r G /C^ ag : 

the matrix element of C mag Oc 


This coefficient is the final result after the use of equations of motion. We prefer to give the coefficient 
in front of the renormalizaton group invariant combination that enters the HQET Lagrangian. This 
combination also determines the mass splitting in the ground state multiplets due to spin orientation. 

Thus, the final coefficient of the matrix element of the chromo-magnetic operator with account 
of equation of motion after taking hadronic matrix elements reads 



-CJr) + 


Ch(r) 

Cmag (/^) 


(40) 


This is a coefficient in front of the matrix element of the renormalization invariant combination 

c ma , g (n)o G (n)- 

In Fig. [6] we plot the mass dependence of this final coefficient. 
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Figure 6: The mass dependence of the coefficient C/^(r) at LO and NLO in the pole mass scheme. 
Color blind: C A = 3, Cp = 4/3. Left panel - the whole phase space, right panel - zoomed image of 
the small mass region 0 < r < 0.1 


Writing again the decomposition of the whole coefficient in a s order 


C-M = Clf + ^ \ C A C, 2 (r) NLO ’ A + C F C, 2 (r) NLO ’ F 


V'G ^ M g p f A ^G ' 

we obtain at the leading order the well known result 

iLO 


PG ' 


(41) 


C -2 = — 3 + 8r — 24 r 2 — 12 r 2 ln(r) + 24 r A — 5r 4 . 


(42) 


The whole expressions are given in Appendix [Cj Here we present the new result at NLO as a small 
r expansion only 


qNLO,A 


2 tt 2 r 


P-g 


17 r 


87 \ 2 Jr 2 7T 2 31 

_l_ r j Q g ^ _ 25 r log(r) —— + — , 


C N 2 LOF = -47r 2 r-19r + 

/4 


3 0 w ow 9 18 

327r 2 y / r 2 57 t 2 91 

---4 r log (r) + 68r log(r) - — - — . 


(43) 


The color blind expansion for QCD (C A = 3, Cp = 4/3) reads 

567 r\/r 

H-— 

9 


1NLO 

f 94 

Pg 

\27 _ 


1 , 


+y(- 


1 2 

H-r 2 


81 


280 

~W‘ 


(44) 


The very large contribution of the -y/r term leads to a very fast change of the coefficient C^ 2 LO from 
its massless limit value with an increase of the charm quark mass. Numerically one finds 

= —3.46 + QIAly/r + r (—2.3 ln 2 (r) + 15.7ln(r) - 109.2) + O (r 3/2 ) . (45) 
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In the massless limit the new result is 


% = 


3+ tM is-9 


31 


c ,'91 Stt 2 
° F 1 72 + 18" 


Note that the Cf part of this coefficient differs from the result given in ref. US]. The difference is 


given by 2 Cq LO and it emerged because 


m 


only the leading order of C 0 coefficient was used for 


subtracting the contribution of the local bfb operator. 


The fi dependence of the prefactor of Oq in Eq. (13) matches the leading order anomalous 
dimension of chromo-magnetic operator [27], such that Cp? is /! independent. 

The end-of-spectrum behavior reads 


qNLO 


Pl (r~> 1) = 40(1-7-) + 


—Cf - C A 


(1 — r ) 5 + 0((1 — r ) 6 ) 


(46) 


for NLO and 

C$(r -► 1) = -4(1 - r) 4 + 0((1 - rf) (47) 

for the leading order contribution. 

The mass parameter of the heavy quark mb is chosen to be the pole mass which is a proper 
formal parameter for perturbative computations in HQET (see discussion in [20]). After having 
obtained the results of perturbation theory computation for the coefficients of HQE, one is free to 
change this parameter to any other [42] . 


6 Discussion of the results 

6.1 The total width 

The radiative corrections are of reasonable magnitude and are well under control for the numerical 
values of the coupling constant for fi ~ 2 — 4 GeV (for the numerical value see, e.g. |43j. This 
provides a clean application of the results to phenomenology. The final quark mass dependence is 
remarkable. It is very fast for small m c therefore the decays into light quarks u for bottom mesons 
and d for charmed mesons should be treated with care. 

The coefficients of HQE have been also calculated in ref. [H] where the analytical computation 
has been performed for the hadronic tensor and the final integration over the phase space has been 
done numerically. Such a setup has advantages for direct comparison with experimental data since 
the experimental cuts in the phase space can be readily introduced. 


19 





We can make a literal comparison with the results of [43J for the total width. Onr result in the 
format of ref. [33] is 

( 2 2 \ 

(1 - 1.7776—)(1 - - (1.9449 + 2.4235 —)^f ) (48) 

7T Z777/^ 7T 777/^ J 

for r = 0.0625 that literally coincides with the results of ref. [33] , 

For phenomenological applications and comparison with experiment within our approach one 
can compute moments of the differential distribution (see, e.g. BSD- It is straightforward to compute 
almost any moment in the invariant lepton pair mass, lepton pair energy or invariant mass of the 
hadronic system. We present few such moments below. 


6.2 Moments of differential distribution 


Note that our computation is organized such that it allows for computation of certain moments of 
differential distribution. We can build up moments over the leptonic pair invariant mass squared 
q 2 , (q = Pi + p u ) and the partonic invariant mass squared (p — q) 2 , p is the momentum of the 
bottom quark and p = m^v. It is possible because we have the leptonic part and the partonic parts 
separately in an intermediate representation of computed diagrams - one can compute the moments 
in q 2 or/and in (p — q) 2 . The total lepton energy moments (the moments in the variable pq) are 
just the linear combinations of those two sets. We present the analytical results for few moments 
at small r expansion for brevity. The analytical expression for the total width is given for further 
comparison with the moments. It reads 


r/r„ = l-8r + CV^| (y-27T 2 -8r(61n(r) + 17) ] + 


^(_3 + 8r+ ^x 

2 mr l 47T 


c A 


27r 2 r 


17r 


+Cp ( —47r 2 r — 19r + 


8n 2 \/r 2 . 7T 2 31 

-b r hr (r) — 25r ln(r)-1- 

3 V 7 K J 9 18 

327 r 2 \/r n, . 57 t 2 91 

-4r In (r) + 68r ln(r)- 

3 18 72 


(49) 


The normalized q 2 moments of the total width with C ^ coefficient are given below. For con¬ 
venience they are normalized to unity at leading order of power, small mass, and perturbative 
expansions. The normalization can be obtained independently. Indeed, the x = q 2 /ml distribution 
in massless limit at LO is given by 


1 dr , 

r - = 2(1 - x)il+2x) . 


(50) 
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The normalization factors for the moments n — 1 — 3 are then N(M%) = {3/10,2/15,1/14}. For 
example, 

l 

N(Mf) = / 2(1 - x)(l + 2 x)xdx = 3/10 . (51) 


First moment (g'/m ^) 1 is 


Ota 


Mf = 1 — 15r + C F ^~ ( 13 - 2vr 2 - r ( 90ln(r) + 42T- + 355 


107r 2 


+ 


he 
2 m 2 


-f + -1 C A 

3 7T 


807b 


25tt 2 260 \ ^ /320tt 2 

^ + 77 ) +Cf {-^ 


+ 


Second moment (g 2 /m 2 ) 2 is 


Ml = 1 - 24r + C F ^ - 2 tt 2 + r (-144 ln(r) - 2 tt 2 - 6813 ^ 


47T 


10 


) 


l l G 
2 m 2 


a. 


-15 + — Ca 

7 r 


Third moment (g 2 /m |) 3 is 


-20tt 2 


a, /1243 


17tt 2 541 


12 


+ 


40 




657T 2 

54 " 


+ 


041 \ ^ _ 77T 2 

—— ) + Cp ( 8 O 7 r v/r H—— 


763 

"36" 


41 


Mn —l — 35 r + C F I - 
3 47T V 


/4 

2m? 


-23 + — ( C A 
7 r 


90 

1127r\/r 


— 27t 2 + r ( —210 log(r) — 


147T 2 


20195 V 

) l + 


18 


+cv 


/ 4487b 


+ 


3 

677T 2 




35tt 2 27217\ 

18” + 1620 ) 
1088429\ 

16200 ) 


(52) 


3 18 

The q 2 moments are very stable and hardly change with n besides the total normalization. Usu¬ 
ally one argues that radiative corrections should increase or decreases depending on the momentum 
flow through the diagram - we see no simple explanation for the change of radiative corrections. 
The moments in partonic variable (p — q) 2 — m 2 are defined through the relation 

((p — q) 2 — m 2 ) n dr 


M h = 

n 


~.2n 

l b 


(53) 


mr r 0 

and have been considered in m- They are given below for n = 1 — 3 analytically within small r 
expansion. 

First moment ({p — q) 2 — m 2 ) 1 is 


M i = C F — ( + h r log(r) + 422- ) + 


cr o 


71r 


91 \ 


he 

2m? 


7T 

7 r 


24 ' 2' 
611r 


+Cv(- 


108 

737r 2 r 457r 

' 108" 


6ooy 

22 29 \ 

—rlog(r)- 




9 


180 ) 

2 


36 


67 , . 7T 2 77 \ \ 3r 1 

-rlog(r)- T + - - 


(54) 
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Second moment ({p — q) 2 — m 2 ) 2 


Mo = C F — | 


0t« 


7T V 432 600 


1% a s 


137 r 


2 ml 7 r 


163 


Ca | T ( —— + 

1 1 18 1080 


+ 72 1 +Cf 


25 In r 3703 \ 


18 


+ 


1080 


J 


+ 


347 \ 
3600J 


(55) 


Third moment ((p — q) 2 — m^) 3 is 



a, / 377 
F 7i V176400 


1^l) + ^( Ca J^ + c f ^L). 

3600 ) 2 m 2 Ti V 16200 1587600 ) 


This set is such that moments vanish at leading order. Therefore one cannot discuss the relative 
magnitude of radiative corrections. Our results for n = 1 — 2 coincide with those of ref. 0®. 

We also compute the relevant moments numerically with full mass dependence for a typical value 
of the mass ratio. For the q 2 moments up to third order we obtain with r = 0.0625 

— 2 

M? = (1 - 1.7776 —)- 3.8898(1 - 0.9206—)-^, 

7r 7r 2 mi 

— 2 

Ml = (1 - 1.6500 —)- 8.9901(1 - 0.6834—)^%, 

7T 7T 2 ml 

— 2 

Ml = (1 - 1.5575—)- 14.394(1 -0.5578—)-^, 

2 7T V 7T 2 ml 

— 2 

Ml = (1 - 1.4847—)- 19.997(1 - 0.4666—)-^. (56) 

6 7T V 7T 2 ml 


Numerically for partonic moments in H = (p — q) 2 — m 2 c with r = 0.0625 one obtains 


— 2 

= 0.0569—+ 0.397(1 - 2.304—)^%, 

7T 7r 2m^ 

Mf = 0.00575— + 0.0554—, 

Ti Ti 2mg 

= 0.00114— + 0.00694—, 

d 7T 7T 2m? 


(57) 


where = C mSig (p)p 2 G (p). 

It is also possible to compute the moments of the lepton energy spectrum that is of interest 
from the experimental point of view. However, here a few more technical problems arise. On the 
one hand the whole set up of the analytical calculation has to be modified, since lcptonic tensor 
has to be taken as a differential distribution rather than fully integrated over the lepton phase 
space. On the other hand there is the question of how to deal with y 5 in dimensional regularization. 
For the cases we discussed here we always have a situation when there is an even (in fact two) 
number of 75 -matrices within the trace over Dirac matrices both in leptonic and hadronic parts, so 
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we simply and consistently use anticommuting 75 . However, in the calculation of the moments of 
the charged-lepton energy one has also consider an odd number of 75 -matrices in the traces, which 
causes an additional complication of the calculation. Nevertheless, with the technology developed 
here, these problems can be tackled and we plan to present a calculation of lepton-energy moments 
in a separate publication. 

6.3 Phenomelogical outlook 

This paper has been devoted to the description of the technical aspects of the calculation of the 
perturbative QCD corrections for subleading powers in the 1/m expansion. Aside from more the¬ 
oretical consideration, such as the discussion of the mass dependence of the various terms of the 
heavy quark expansion, such a calculation has a variety of phenomenological applications, of which 
the most prominent one is its application to inclusive scmileptonic b —> c transitions. 

These decays are currently believed to be the most precise method to determine the CKM matrix 
element V c b- In this method, V c b is extracted form the heavy quark expansion for the total rate, 
while the heavy quark expansion parameters /i n , pr> etc. are extracted from the moments of the 
differential rates. Based on this methodology, the theoretical uncertainty in V c b has been reduced to 
a level below 1 %, while the total uncertainty (including the experimental as well as the uncertainty 
in the extraction of the heavy quark expansion parameters) is at the level of 2%. The current 
extractions of V c b do not yet include the a s /j, G contributions, which are parametrically the largest 
missing pieces in the analysis. 

From the experimental side, the lepton energies cannot be measured to arbitrarily low values. 
Thus either an extrapolation is necessary or one has to include a cut into the theoretical predictions. 
Since an extrapolation involves a model dependence, it is more favorable to include a lepton-energy 
cut into the theoretical prediction. 

However, unlike in the numerical study of (26) M], such a cut cannot be implemented in an 
analytical calculation, at least not exactly. Thus in order to make phenomenological use of the 
analytical calculation one needs to take into account effects of such a cut, which needs further 
study. We plan to return to this in a separate publication. 
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Appendices 

A Master integrals 


Here we present the results for master integrals entering our calculation in dimensional regularization 
with D = 4 — 2e. The general integral 

S(a,b;m 2 ,p 2 ) = J (m 2 _ p)a(_( p _ k yy ( 58 ) 

develops a cut at p 2 = s > m 2 with a discontinuity p(a, 5; m 2 , s) that is 

p(a, b ; m 2 , s ) = -— T (S(a, b; m 2 , s + fO) — S(a, 6; m 2 , s — iO)) . (59) 

2m 

It is a spectrum of a general sunset diagram [30] 

, , 2 , T(D/2-b) 

• S) = r(a)r(b)r(D - a - 2b + 1) 

z D-a-2b s Dl2-a-b 2Fi ^^ ^ X £) _ a _ + 1; «) (60) 

with z — 1 — m 2 /s. Here 2 Ti(a, 6 ; c; 2 ) is a hypergeometric function. I 11 our case m = m c and 
s = m 2 . 


A.l Master integrals at LO: two loop 

At LO there are two master integrals. In both cases it is a two loop sunset with one heavy (m c ) 
line. The internal massive line can be a normal one or doubled which is denoted by a dot on it, see 
Fig.0 



Figure 7: Two-loop master integral. A dotted line indicates one additional power of the propagator. 


The closed form for a master integral with a normal line is 


M 00 = 'S'(l) 1; 0, —l)p(l, 2 — D/2;m 2 c ,m 2 b ) 


(61) 
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where ^(l, 1; 0, —1) is a scalar massless loop that is expressible through h-functions 

We usually set mj = 1 in the computation. The e-expansion of this integral can be obtained with 
the program HypExp or independently. At the leading order of e-expansion one has 


I 777^ 

Mm = - + m 2 In (m 2 ) - — + OU i 


(63) 


with m = m c and mb = 1. 

The second master integral (dotted) belongs to the same class of sunsets and can be obtained 
as a derivative in m c 

Moi = ~d^ M °°- ( 64 ) 

The closed form for this dotted leading order master integral is 


Afoi = S (1,1; 0, -1) p( 2,2 - D/2; m 2 , m 4 2 ). 


(65) 


A.2 Master integrals at NLO: three loop 

At NLO there are master integrals that are factorizable, of sunset-type, and nontrivial. 

A.2.1 Factorizable integrals 

The factorized master integrals contain a closed massive loop that can be either of charmed quark 
or bottom quark. 



Figure 8: Factorizable three-loop master integrals. 


These master integrals are: 


with T 0 (m) being a massive tadpole 


Mi i — To(m c )Moo 


T 0 (m) = m D ~ 2 Y{l - D/2) 


( 66 ) 


(67) 
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and the other one 


( 68 ) 


M\2 — T 0 (m c )M 0 i. 

The master integrals with a 6-quark tadpole are 

M 4 i = T 0 (mb)M 00 

and the other one 

M 42 = T(m b )M 0 1. 

In actual computation the bottom quark mass is set to unity. 

A.2.2 Sunset-type integrals 

Non-factorizable but still simple master integrals M 2 i,M 2 2 are of the sunset type. 


(69) 


(70) 



Figure 9: Sunset-type three-loop master integrals. 


The normal one is given by the basic integral 


A /21 = 5(1,1; 0, -1)5(1,2 - D/ 2 : 0, —l)p(l,3 — D,m*, mi) 


( 71 ) 


The dotted one is its derivative in loop (charmed quark) mass 


d 

M22n — —^dT 2 i 

am* 


(72) 


which is again a three loop sunset 


A/ 21 = S(l, 1; 0, —1)5(1,2 — D/2;0,-l)p(2,3 - D, m/. m/l . 


(73) 


A.2.3 Nontrivial master integrals 

There two nontrivial master integrals that can be chosen in a varity of ways. We define the first 
nontrivial master integral N p as a sum of left (dotted at bottom line) and right (dotted on charm 
line) diagrams in Fig. ??. In words, this can be expressed as N p = dot.m b + dot.m c . 
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Figure 10: Nontrivial three-loop masters. 


We managed to compute the e-expansion of N p up to the necessary order. It reads 


N v = NZ° + eN, 


tNLO 


(74) 


with 


Np = -2(1 - r) - (1 + r) log(r) 


(75) 


and 


N. 


NLO 


4a Jr ^4(Li 2 (— yfr) - Li 2 (\/r)) + 7r 2 + 2 ln(r) In 
+ Q - 0 ln 2 (r) - 3 (r + 1) ln(r) 

+4(r + 1) ln(l — r) ln(r) + 8(1 — r) ln(l — r) + 14(r — 1). 


(76) 


These are master integrals entering partonic contribution for the total width and C v coefficient. 

At NLO one more master integral appears to be necessary for the Cq coefficient. It is represented 
by the difference of left and right diagrams in Fig. ??. In words, this can be expressed as N p = 
dot.rrib — dot.m c . We need only the leading term of its e-expansion that reads 


27T 2 

N m = (1 — r) ( —4Li 2 (r) H——-41n(l — r) ln(r) + 21n(r) ) — 2rln 2 (r) 


(77) 


A.3 Master integrals in massless case 

We have calculated all quantities in the massless limit independently. The reduction procedure and 
master integrals have been obtained independently as well. In massless case master integrals can 
be found in a concise form. These master integrals are represented by Feynman diagrams given in 


Fig. 11 


At leading order there is one master integral 


M 0 ° 0 = 5(1,1; 0, -1)5(1, 2 - D/2; 0, -1) 


sin( 27 T£) 


7r 


(78) 


At NLO in massless case there are three master integrals: 
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Figure 11: Massless two- and three-loop master integrals. 


a) factorizable integral 


M° n =T 0 (m b )M° M - 


(79) 


b) sunset integral 

= 5(1,1; 0, -1)5(1, 2 - D/2 ; 0, -1)5(1, 3-D;0, -1) ; (80) 

7r 

c) complicated integral 

P(T _ c-i 2 

N a = —5(1,1; 0, ~ 1 ) r(2 b r(3 _ 3e ) ^(fe 1 - £ .!}; (3 - 3e, 2 - e}, 1). (81) 


B ( coefficient at NLO with full mass dependence 

The expression for the coefficient C v is 

C N L ° = ( 3 Li 2 (r) - ^tt 2 ) (1 - 16r 2 - 3r 4 ) - ^-(1 - r) (25 - lOllr - 1487r 2 + 189r 3 ) 

+—r (12 + 450r + 4r 2 + 45r 3 ) ln(r)-(1 — r) (ll + Hr + 83r 2 — 45r 3 ) ln(l — r) 

6 6 

3 

+ -r 2 (4 + r 2 ) ln 2 (r) + 2 (l — 30r 2 — 3r 4 ) ln(l — r) ln(r) 

_)_ 8 r 3 / 2 (i _)_ 3 r ) ^4LiJ — 7 T 2 — 2 In ^ ln(r)^ (82) 

where Li^ = Ll 2 (\/r) — Li 2 (— y/r). 
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C Cq coefficient at NLO with full mass dependence 


Here we give results for the C^ 2 LO coefficient. At NLO we give both color structures. 
The Ca color structure coefficient of a^/V reads 


(jNLO.C'a 


a = —(1 - r) (156 - 4081r - 354r 2 - 405r 3 ) 
Pg 108 v ’ 

H—( 6 Li 2 (r) — 7T 2 ) (l — 6 r + 24r 2 — llr 3 ) 

9 


(1 


54r 


(15 + 20r - 196r 2 - 292r 3 - 27r 4 ) ln(l - r) 


1 2 

—r (786 + 972r + 131r 2 — 27r 3 ) ln(r) — - (l + 9r — 93r 2 + 19r 3 ) ln(l — r) ln(r) 


54 

H—r (9 — 33r + 5r 2 ) ln(r)^ 
9 


+ !^ 1/2 (l “ y r ) ( 4Li 2 - vr 2 - 2 ln(r) In Q ^ ^ 


The Cf color structure is 


(83) 


(j NLOfip 


P-G 


= - —(t - r) (321 - 13747r + 5421r 2 - 3807r 3 ) 

+ — ( 6 Li 2 (r) - 7 T 2 ) (5 + 72r - 72r 2 - 88 r 3 + 45r 4 ) 
18 


(! ~ r ) 
54 r 


(12 - 19r + 917r 2 - 1795r 3 + 585r 4 ) ln(l - r) 


+ — r (1500 - 330r + 2668r 2 - 585r 3 ) ln(r) 

2 

H— (ll + 54r — 48r 2 — 94r 3 + 45r 4 ) ln(l — r) ln(r) 
9 

-—r (72 + 60r — 112r 2 + 45r 3 ) ln(r) 2 

18 

+ y (! - ^ r)r x/ 2 ^4Li2 - vr 2 - 21n(r) In Q ^ 


(84) 
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